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We calculate the first-order perturbative correction to the transition temperature T c in a 
superconductor with both non-magnetic and magnetic impurities. We do this by first evaluating 
£Nj ' the correction to the effective potential, f2(A), and then obtain the first-order correction to order 

parameter, A, by finding the minimum of Q(A). Setting A = finally enables T c to be evaluated. 
T c is now a function of both the resistance per square, Ptu, a measure of the non-magnetic disorder, 
and the spin-flip scattering rate, l/r s , a measure of magnetic disorder. We find that the effective 
^ ■ pair-breaking rate per magnetic impurity is virtually independent of the resistance per square of the 

film, in agreement with an experiment of Chervenak and Valles. This conclusion is supported both 
by the perturbative calculation, and by a non-perturbative re-summation technique. 
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I. INTRODUCTION 



Many experiments^! performed on homogeneous disordered thin film superconductors have shown that super- 
i-^j | conductivity is suppressed by increasing disorder, as measured by the normal state resistance per square, Rn. The 
£^ ■ majority of the data is of the transition temperature as a, function of the resistance per squarje 1 -T c (i? n )B~B, although 
there is some data on the upper critical field, H c2 (T, i?n)&EjQ, and the order parameter A(R a p , u. The main data to 
be explained thus consists of T c (Ra) curves for different materials. To see how disorder might affect T c , consider the 
mean field equation, 



T c a = 1.13(xi_Dexp 



i 



iV(0)(A-/i* 



(1) 



where A is the attractive BCS interaction mediated by phonons of energy less than the Debye frequency, Wjj, 
| is the Coulomb pseudopotential, the effective strength of the Coulomb repulsion, and N(0) is the single particle 
density of states at the Fermi surface. Obviously disorder could affect T c by changing A, /!*, and iV(0). In particular 
the diffusive motion of electrons caused by the disorder is known to lead to an increased effective strength of the 
"p^ . Coulomb interaction, as the screening is less efficient than with ballistic electrons, and this leads to an increase in 
/i*, and a decrease in j'V(O). Calculating thaJksjt-order perturbative correction caused by the disorder shows that 
we must consider all these processes togetherEElEJ. This is because the disorder-screened Coulomb interaction has a 
low-momentum singularity which leads to the separate effects being large; however, when they are added together 
this singularity is cancelled, and the actual effect is much smaller than might be naively expected. The final result 
^ i has the form 

\T c0 J 3R \2TrT c oTj 7 1 ' 

where Rq — 2-Kh/e 2 sa 162fcf2 and r is the elastic scattering time. We see that this curve is essentially "universal", 
depending on only a single fitting parameter, (3 = In (1/2ttT c qt). Experimentally it is found that T c {Rn) curves from 
a wide variety of materials fit well to this equation, or extensions of it that allow for stronger disorder. The simplest 
extension simply consists of replacing T c q by T c on the right-hand side of Eq. (2), which leads to the cubic equation 

x= t -{{3 + xf, (3) 

where x — In (T c q/T c ) and t = R u /Rq. This equation shows unphysical reentrance at strong disorder, an artefact which, 
is removed by either a renormalization group treatmentli2l, or the use of a non-perturbative resummation techniquetHI 
to yield the formula 
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(4) 



The fact that most data can be fit to a single curve is pleasing in that it shows that the basic ingredients of 
our theory - disorder, BCS attraction and Coulomb repulsion - are correct. However it does not allow analysis 
of the sensitivity of the theory or experiment to changes in details of the system, such as the exact form of the 
phonon-mediated attraction. Moreover there are other theoriesliS which posit the importance of such details which 
give predictions that are equally in agreement with experiment. We see that the X" c (-/?□) curves alone are not enough 
to allow consideration of the relative merits of different theories. What we would like to do is to add some additional 
parameter to the experimental system to give a whole new set of data p-for example a family of T c (Ru) curves for a 
single material as this new parameter is altered. Chervenak and Vallcaij have recently performed an experiment of 
this type in which magnetic Gd impurities are added to thin films of P6q.9-ZKo.i- This introduces the new feature of 
spin- flip scattering to the system, which is measured by the spin- flip scattering rate, l/r s . The task of the theorist 
is to now make predictions for T c as a function of both Ra (the measure of non-magnetic disorder), and l/r s , (the 
measure of magnetic disorder), and to compare these to experiment. 

In this paper we calculate the first-order perturbative correction to the transition temperature, T c , of a super- 
conductor with both non-magnetic and magnetic impurities. The model used consists of a featureless BCS attraction, 
— A, and a Coulomb repulsion, Vc(q), between electrons which scatter off non-magnetic and magnetic impurities. The 
model is the simplest one that contains the essential physics, and its shortcoming of not considering the details of the 
attractive interaction is offset by the fact that we can consider all processes to a given order of perturbation theory. 
This is an important consideration in view of the cancellation of low-momentum singularities in the screened Coulomb 
potential discussed in the opening paragraph. In fact, an obvious question is whether this cancellation persists in the 
presence of magnetic impurities. We find that this is indeed the case, and so the details of the screened Coulomb 
interaction are removed, leading to a "universal" form for T c (Ru, 1/t s ). 

The main result of the paper is that the pair- breaking rate per magnetic impurity, a'(i?n), defined by 



is roughly independent of Ra except near the superconductor-insulator transition, in agreement with experiment. This 
is confirmed both by first-order perturbation theory, and also by a non-perturbative resummation technique which we 
introduce to remove concerns about reentrance problems at stronger disorder. This agreement of the two theoretical 
approaches with each other and the experimental data gives us confidence in our results. ■— . 

To calculate the correction to T c we use a collective mode formalism derived in a previous papeiEd (which we 
refer to as I from now on) on the suppression of A by non-magnetic disorder. The introduction of magnetic impurities 
means that we have to modify the formalism somewhat, and so we include most of the derivation in this paper. The 
method used in this paper to evaluate the correction to T c proceeds in three stages. First we find the first-order 
correction to the grand canonical potential, f2i(A), of the superconductor due to fluctuations of its collective modes. 
Then by minimizing the total grand canonical potential, r2o(A) + f2x(A) with respect to the order parameter, A, we 
obtain the first-order correction to the order parameter self-consistency equation. Finally by setting A = we obtain 
the first-order correction to transition temperature T c . The method we use has the advantage that it is impossible to 
"miss diagrams" since there is only one diagram in the Qi(A) calculation, and we also obtain the equation for A at no 
extra cost. The equations for T c and A must reduce to those of I when we set spin-flip scattering to zero, providing 
a useful consistency check. A key result of the calculation in I was that the singularity in the screened Coulomb 
potential persists below T c , and that this singularity is cancelled in the formula for the suppression of A in a similar 
manner to the cancellation in the formula for T c . Therefore an important question is whether this singularity and 
cancellation remains when magnetic impurities are added. We show that this is indeed the case, and moreover that the 
cancellation is due to gauge invariance, and will occur in first-order perturbation theory in the presence of any kind of 
impurity scattering. In other words, it is not possible to obtain stronger suppression of the transition temperature by 
introducing some exotic scattering mechanism. It is also reassuring to know that an otherwise mysterious cancellation 
between diagrams has its physical origin in gauge invariance, and we hope that similar arguments may be applied to 
show that the result holds to all orders in perturbation theory. 

The outline of the rest of the paper is as follows. In section II we derive the matrix formalism for superconductors 
with magnetic impurities, and the collective mode approach we will use. We derive the RPA screened bosonic 
propagators, and show that low-momentum singularities persist in the screened Coulomb propagator below T c . In 
section III we derive the first order perturbative correction to the grand canonical potential Qi(A), and from this the 
correction to the order parameter A. In section IV we set A = to obtain the correction to transition temperature T c . 
In section V we calculate T c numerically using both the perturbative results of section IV and a recently developed 
non-pertubative technique, and compare to experiment. 
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II. DERIVATION OF THE 4x4 MATRIX FORMALISM 



Superconductivity with Magnetic Impurities: 

We consider a system of electrons that scatter off static non-magnetic and magnetic impurities, and interact with 
each other via the long-range Coulomb interaction and the BCS attraction. The scattering from static impurities is 
described by Hamiltonian 



ff e _< = W dx^+(x) J 

a/3 J [ 



2m 



(6) 



where ip£, ip a are the electron creation and annhilation operators, uo(x — Xj) is the impurity potential at x due to a 
non-magnetic impurity at x^, Sj is a magnetic impurity spin moment at xj, and J(x) is the electron- impurity exchange 
coupling. The potential and spin-flip scattering rates are then given by 



TQ 
1 



2TTN(0) ni \u \ 2 
27rN(0) nj J 2 S{S + i), 



(7) 



where uo is the Fourier transform of ito(x), and J is the Fourier transform of J(x), both assumed independent of 
momentum, is the non-magnetic impurity density, and rij the magnetic impurity density. 
The Coulomb repulsion between electrons is described by Hamiltonian 



H C = Y, / dx / dx>+(x)^(x)^-^+(x')#(x') 

a/3 |X X 1 



(8) 



leading to a bare Coulomb propagator that is just the Fourier transform of the above potential. 
The BCS attraction is described by the Hamiltonian 



Hbcs = / dx ^( x )' ( /'a(x)^(x)'0/3(x), 

a/3 J 



(9) 



corresponding to an instantaneous contact interaction — A<5(x — x'). 

Having introduced the model—Hamiltonian we need to describe the system, we discuss the standard four- 
dimensional matrix representation!!^ needed to describe a superconductor with magnetic impurities. We need four 
components to describe the two spin degrees of freedom, and the two types of correlation - the usual particle-hole 
correlation < ipiip^ >> amd the anomalous pairing correlation < ipfip± >■ We introduce the four-dimensional vector 
operator 



* = 



1>t 



"01 W 



(10) 



with matrix propagator 



< >= 



< > 



< W£ > 

< ip+ip+ > 

< ip+ip+ > 

In the normal state the temperature Green function is 



< > 



< i>ii>i > < i>if \ > 

< Vj. > < V'fV'T > 



(11) 



(12) 



where z = iuj, to = (2n + 1)ttT is a Fermi Matsubara frequency, and the T; and cr,; are Pauli matrices operating on 
different spaces. The Oi operate in the usual spin space, whilst the Tj operate in the Nambu (electron-hole) space. The 
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diagrammatic rules are then the same as in the normal state except for the matrix structure of the Green function, 
and the presence of matrix t 3 cq at each interaction or impurity vertex due to the electron density operator being 
written in the form 

p = V^Vt + V'+V'l = i* + T3(To*. (13) 

The pairing correlations in the clean superconductor can then be taken into account self-consistently as shown in 
Fig. (la). Making the ansatz £ = AT1CT3 for the self-energy, the Green function for the pure superconductor becomes 

r , (u x 1 z + e k T 3 a + Atict 3 

G 0(fc, Z) = -r = g ^2 ' ( 14 ) 

z - s k T 3 a - Ati<7 3 z 2 -e%-A 2 
and the diagram of Fig. (la) gives self-energy 

S = -\TY N(0) [ de k T3ao{z ± £ 2 kT3ao 9 + 

^ J e%- z 2 + A 2 

= iV(0)AAr 1( r3T^ 1 (15) 
~f + A 2 

which gives us the usual BCS self-consistency equation 

We can treat the presence of non-magnetic and magnetic impurities by including an extra self-energy diagram 
to describe the dressing of the electron line by impurities as shown in Fig. (lb). We then make the ansatz that the 
pairing energy has the form S p = Ati<j 3) and the impurity self-energy has the form ^ imp = — (z — z) + (A — A)ti<73, 
so that the Green function for the dirty superconductor is 

Go(k, z) = — — 2 , (17) 

z 2 - e\ - A 

which is just the Green function for the clean superconductor with z, A, replaced by z, A, respectively. Since the 
impurity line has the form 

r ° = r, ,! n ^ T 3O0 ® T 3 <T + 1 [ToCTl <E) T(](7l + To<7 2 <E) Tq<T 2 + T 3 (T 3 <g> T 3 <T 3 ] , (18) 

we obtain the self-consistency equation for ~z = iZJ and A, 

The diagrammatic definition of the pairing energy, E p , leads to the same self-consistency equation for A as in the 
pure case except that u), A, are replaced by 57, A. In the absence of magnetic impurities - i.e. l/r s = - we see that 
lu/A = uj/A, and the equation for A is unchanged. This is Anderson's theoremlij that superconductivity is unaffected 
by non-magnetic impurities at mean-field level. In the presence of magnetic impurities we see that ZU/A ^ uj/ A, and 
if we define u =ZJ/ A, ( — 1/t s A, the problem reduces to solving the equation 

" =u fi_ c _^ = y (20) 

A v W/ 

The self-consistency equation for A then takes the form 

1 = 7V(0)AT^^=, ( 21 ) 
and in particular if we set A = we get for T c , 
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1 = N(0)XT c 



l/r s 



(22) 



Subtracting c£L the equation for T c q, the transition temperature in the absence of magnetic impurities leads to the 
famous resultO 



log 



1 



2ttT c t s 



(23) 



Collective Mode Formalism and RPA: 



The idea of the collective mode formalism is to treat the screened interactions in the system as bosonic collective 
modes. The relevant bosonic operators are order parameter amplitude and phase, and electron density. These are the 
only modes that are coupled by a bare interaction - the BCS interaction for order parameter amplitude and phase, 
the Coulomb interaction for electron density. The main advantage of this approach is that we are able to treat order 
parameter fluctuations and the Coulomb interaction on an identical footing. This procedure may be formally carried 
out within the path integral theory of superconductivity by decoupling the two four-fermion interaction terms with 
the introduction of appropriate collective variables. This is discussed in detail in the paper of Eckern and Pelzerlia. 
The end result is that there are three effective bosonic modes, order parameter amplitude and phase and electronic 
density, and each can be written in the form 



O, 



2 



(24) 



with matrices Ma = t\Uz for order parameter amplitude, = T2CT3 for order parameter phase, and M p = T3CT0 for 
electronic density. Interactions occur by exchange of these collective modes, and so the effective interaction potential 
is now a 3 x 3 matrix. The screened interaction is found from the equation 



£ 

kl 



(25) 



V° = 




as shown in Fig. (2). Here ITy is the polarization operator and VR is the bare interaction matrix which is given by 

\ 

-A/2 , (26) 
Ane 2 /q 2 ) 

the BCS attraction being split equally between the two order parameter modes. The only new diagrammatic feature is 
that any of the matrices Mi can now appear at an interaction vertex, corresponding to interaction with the collective 
variable described by that matrix. 

In order to carry out the calculations in section III, we need the impurity dressed RPA polarization bubbles, Ily, 
shown in Fig. (3). To evaluate the polarization bubble 11^ we must first evaluate the geometric series 

n = s + sr s + sr a sr s + . . . , (27) 

where T is the impurity line defined in Eqn. 

S 



|18[ ) , and S is the momentum sum of a direct product of Green functions 
q,ico + iCl). (28) 



5^G(Jb,tw) ® G{k 

k 



Since we do not need the complete matrix structure of II, but just its traces with two matrices from the set Ti<T3, 
T3C0, we actually evaluate the impurity dressed vertices IX, which have one matrix from the above set inserted 
between the two terms of the direct product in II. These satisfy the equations 



n, 



Si + ST IL 



(29) 



where Sj is obtained by inserting the matrix Mj between the two terms of the direct product in S. These equations 
are solved in Appendix A to obtain the results 



n 



A 



tt7V(0) 
D+ 

ttN(0) 
D~ 
nN(0) 

D- 



UU' 



uu — 1 



UU' 



UU' 



■ uu 



UU' 
i(v! — i 



UU' 



UU' 
■ uu' - 



UU' 



-T10-3 



T103 



UU' 



1 i(u' 
- + 



UU' 



-T1CT3 



(30) 
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where U = \Ju 2 + 1, u' = u(lu + fi) = u(o/), {/' = \/u 12 + 1 and 



ZV + AJ7 + At/' 



1 /" uu' =F 1 
UU' 



(31) 



We can finally obtain the non-zero polarization bubbles by inserting the second matrix from the set ti<T3, T2CT3, 
T3(To into Ilj, taking the trace, and recalling the factor —1 for a fcrmion loop. This yields 



U AA (q,n)^TrN(0)Tj2 
n 00 ( g ,Q)=7riV(O)T^; 



UU' + uu' - 1 



£/[/' + W + 1 



Dq 2 + AU + AU' 



UU' - uu' + 1 



UU' 



iW ? ,fi) 



-tt7V(0)T^ 

-7dV(0)T^ 



UU' 
UU' -uu' -1 



Dq 2 + AU + AU' - 



UU' - uu' - 1 



UU' 



UU' 



UU' 



Dq 2 + AU + AU' - 



1 



UU' -uu' -1 
UU' 



Dq 2 + AU + AU' 



UU' -uu' -I 



UU' 



U p ^(q,Q). 



(32) 



We note that if we set l/r s = in Eqn. (|32|) we will obtain exactly the results found in I, as of course we must. The 
screened potentials Vij are then given by 



V 



where 



(-A^ + riAA)- 1 

[(zvciq))- 1 +u pp yv -u^/v 

n^ p /x> (-A-i + n^vuj 



v = (-a- 1 + n 00 )[(2i/ c (g))^ 1 + n pp ] + 



(33) 
(34) 

n„w, is a manifestation 



The coupling between phase and density fluctuations caused by the non-zero value of n^ p 
of gauge invariance. 

We can next show that the propagators V^, V$ p and V pp all have a low-momentum singularity of the form 
l/q d ~ 1 for all non-zero frequencies and all temperatures. In other words, these propagators have the same long-range 
behavior as the unscreened Coulomb potential. An analogous situation occurs for the disorder screened potential 
in the normal metal, where the singularity is known to strongly affect the properties of the system. To show the 
existence of the singularity we simply need to show that the denominator V vanishes at q = for all ft ^ and all 
T. Since Vc(q) ~ l/q d ^ 1 , we need only prove that 



[-A- 1 n 00 (o, o)]n w (o, n) + n 0p (o, ft) 2 = o 

This is proved in Appendix B where we show that 

n 



-A" 1 + 11^(0, fi)] 



2A 



iWo,fi) , n 0p (o,fi) 



iWo.fi), 



(35) 



(36) 



by direct calculation. We also show that this result is guaranteed by gauge invariance and is therefore true no matter 
which scattering mechanisms we include. 



III. FIRST ORDER CORRECTION TO GRAND POTENTIAL AND ORDER-PARAMETER 

SELF-CONSISTENCY EQUATION 

In this section we evaluate the first-order perturbation correction to the grand potential, fi 1 (A). By minimising 
the sum fio(A) + fii(A) with respect to A we obtain the corresponding correction to the order parameter self- 
consistency equation. This method was first used for the system with only non-magnetic impurities by Eckern and 



G 



PelzertS, and we choose to use it as it involves the smallest number of diagrams. The same result for T c can also be 
obtained using the Eliashberg diagrams for A shown in Fig. (5), or the pair propagator diagrams shown in Fig. (6). 

The diagram for the first order correction to the grand potential simply consists of the "string of bubbles" diagram 
shown in Fig. (4) Since the polarization bubbles in this diagram are just those evaluated in the previous section, we 
have all the information we need to derive f2i(A). The only thing we need to remember is the extra symmetry factor 
of 1/n required for the diagram with n bubbles. So whereas previously the RPA equation involved the series 



V = V + V a UV + V UV a UV a 



it now becomes 



1 



1, 



fii = v u + -v^n^n + -v a uv a uv u 



...HV 1 -!!)- 1 , 

= -iog[v 1 -n] + io g [v 1 ]- 



(37) 



(38) 



After summing over all the internal variables - momentum q, Bose Matsubara frequency Q and the three bosonic 
modes - we end up with the final expression for fii, 

fii = E i lo g (~ A_1 + n AA(9, 0)) + log [(-A" 1 + iWg, tymVciq))- 1 + U pp (q, ft)] + Il^q, ft) 2 ] } . (39) 



To proceed we need to minimise the total grand potential, 

d 



— {ft (A) + ft 1 (A)} = 0. 



The mean-field grand potential, fto(A), is given by 

ft (A) = iV(0) A " 



w k 

and after taking the derivative of fto(A) with respect to A, we see that Eqn. (^) takes the form 



A dA ' 



(40) 



(41) 



(42) 



The next step in the procedure is to evaluate <9fti(A)/<9A. From Eqn. ( p9[ ) we see that d/dA will be acting upon the 
Uij to give 



ran 



AA 



\ dA 



V AA 



an, 



OA 



'■Vs. 



' dA 



-V A 



dli 



dA Vpp 



(43) 



From Eqn. (32) we see that the d/dA can act either on the coherence factor or on the denominator in the expression 
for n^-. Acting on the denominator gives a result proportional to the denominator squared, corresponding to the 
two-ladder diagrams of Fig. (5a) in the Eliashberg approach. Similarly acting on the coherence factor leads to the 
one-ladder diagrams of Fig. (5b). We note that the explicit evaluation of the three-ladder diagrams of Fig. (5c) will 
give a zero result. 

The only difficulty in taking the derivatives of the polarization bubbles, n^-, with respect to A is that the quantity 
u(lu) present in all these equations satisfies the transcendental Eqn. (|20|). In Appendix C we evaluate the derivatives 
of the Uij to obtain the results 



dU AA 
dA 

dU^s 
dA 

dn OD 



dA 
dA 




U 3 

u 3 
u 3 
u 3 



1 


u(u' + u) 
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U 3 U' 


D + + 
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u(u' — u) 
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U 3 U> 
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u(u' — u) 
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U 3 U' 


L>_ 
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u{uu' + 1 


) 1 



- 1- 



UU' 



uw 



UU 1 



1 

u 



u 2 



IP 



= 1 



u 2 



u(u' — 


u) 


UU' 




u(u' + 


u) 


UU 1 




u(u' + 


a) 



U 3 U' 



UU 1 



1 UU 

u(u' + u) 
UU' 



1 

Dl 

1 

w 

1 

D 2 
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These formulas together with Eqn. (^3|) lead to our final result for the first order correction to the order parameter 
self-consistency equation: 



1 



N(0)X 



^EtW t E^EE 



Q q 



1 1 u(u' + u) 

a^dZ U 3 U' 
1 1 



1 1 1 
ADlU 



U 



1 - 



u 3 

u(v! — u) 
\- — 

uw 



uu' — 1 
UU' 



A 2 D- 
1 1 



u(u' — u) Tr ^, 2u(uu' + l) Tr . ^, u(u' — u) Tr . „, 
U3U , V <t>4>^ fi ) " Jpfj, - V 4>p(*> n ) + U3U , V pp(i> n ) 



V AA (q,ft) 



u 



u{u' — u) 
UU' 



UU 



1 \ 2(v' — 7/1 



UU' 



v , 



(45) 



The above formula is valid for all temperatures < T < T c , but we are usually interested in the special cases 
T = and A = (i.e. T — T c ). In these two cases the sum over oj on the LHS can be performed analytically to yield 
the two simple forms 



log 



A(0) 
A o (0) 



= ^EE^E- 



n q 



(46) 



n q 



Having noted the presence of the l/q d 1 singularities in the potentials V^, V$ p and V pp , we should now see 



whether the terms in Eqn. (45) containing these singularities cancel out. If we go back to Eqn. (p9[ ) for the correction 
to the grand potential, we see that the term V that goes as q d ~ x is inside a logarithm. Since q d_ ^occurs as a product, 
we can simply take off the term log {q d ^ 1 ) 1 and upon differentiating with respect to A should get zero. In other words 
we naively expect no singular term in Eqn. ( |45| ) . However this is not quite correct since to prove that the denominator 
T> vanished at q — we needed to replace A -1 using the mean-field self-consistency equation, Eqn. (f2l|). We note 
that although the two sides of Eqn. ( |2l|) are numerically equal in the mean-field case, their dependences on A differ 
- A -1 gives zero under d/dA, whilst 1/U does not. This discrepancy then leads to the only singular term in Eqn. 
([45|), which may be written 



In 



A\ 1 1 



1 - 



U 3 



T EE^(^)- 



(47) 



Since this term tends to half the pair propagator contribution to the suppression of T c when we let A — > 0, we interpret 
it as the phase fluctuation contribution. It is singular because of the Mermin-Wagner-Hohenberg theoremt2l which 
tells us that we cannot have broken symmetry states in 2D systems at finite temperature. In the following we will be 
mainly interested in the correction to T c due to Coulomb interaction and so will ignore this term. 



IV. FIRST ORDER CORRECTION TO THE TRANSITION TEMPERATURE 



We can now evaluate the first order correction to the transition temperature by linearizing the order parameter 
self-consistency equation with respect to A. The former can also be obtained directly from the normal state by 
calculating the pair propagator L(q, f2) to first order, and looking for the instability at q = £1 = 0. L is given by 



L- 1 (q,n) = \- 1 + p( q ,n), 



(48) 



where P(q,tt) is the pair polarization bubble. The zeroth order polarization bubble Po(q, O) is shown in Fig. (6b) 
and leads to the mean-field result 



L-\q,n)^N(o) 



N(0) 



log 



T 



T, 



(-00 



1 



Dq 2 



mi 



log (i; 



2itTt s 
1 

2ttTt s 



AttT 

D q 2 + \n\ 



1> 



4nT 



1 



2ttTt s 



(49) 
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where T c qq is the BCS transition temperature (the mean field value in the absence of magnetic impurities), and T c o 
is the mean field value for the system with magnetic impurities. A correction to the polarization operator SP(0, 0) 
will lead to a change in the transition temperature, which is defined as the temperature at which the denominator of 
L becomes zero, given by 



T c \ _ &P(0,0) 
N(0) 



J-cO 



(50) 



If we look at Fig. (6) we see that there are 7 diagrams which contribute to the first order correction to T c . We will 
set A = in the order parameter result of Eqn. ( |4q ) to get the transition temperature equation, and we will be able to 
identify the contribution that comes from each of the Pi diagrams. When we set A — ► 0, then Am — > (|w| + l/r s )sgn(w); 



AU \lu\ + l/r s ; V pp 2V C ; V AA and V+ 



-L; V<h p — > 2H c f >p LVc- The coherence factors then become Heaviside 



functions that set the relative signs of the frequencies 

, uu' ± 1 



UU' 
uu'±l 
UU 1 



26(-u(u + Cl)) 
29(lo(lu + Cl)). 



(51) 



The two denominators, D±, both become Dq 2 + \Cl\ for w, w + £1 of opposite sign; Dq 2 + \2w + Cl\ + 2/r s for w, + 
of the same sign. Making all these substitutions leads to 

1 1 2 



P 1 = -t:N(0)tJ2tJ2J2 



SI q 



1 



P 2 =nN(0)TY,Tj2J2 



(\cu\ + l/T s ) 2 (Dq 2 + \Cl\) + 1/Ts) (Dq 2 + |fi|)2 

1 1 -VcfofiMwCw + n)) 



Vfc(g,fi)«(-w(w + n)) 



n q 



(M + l/r s ) 2 Dg 2 + \2uj + Cl\ + 2/t. 



F 3 = -7riV(0)T^r^^ 
P 4 = -7TiV(0)T^T^^ 

p 5 = -7riV(o)r£r££ 

P 6 =7riV(0)T^T^]T 



1 



(H + i/T a )(|w + n| + i/T.) 
i 



2/t s 



(Dq 2 



mi 



vb(g,n)e(-w(w + n)) 



(|w| + 1/t.)(|w + fi| + 1/t s ) (Dq 2 + \2u; + Cl\ + 2/r s ) 



(Dq 2 + \Cl\) 2 _ 

V c (q,Cl)9(u(uj + Cl)) 



1 



(M + i/t s ) 2 
i 



i 



M - l/r s 



(H + l/ra^W + ini) 



(Dq 2 + \2iu + Cl\ + 2/t.) (Dq 2 + \2lo + Cl\ + 2/t s ) 2 
L(q,Cl)6(-cj(u + Cl)) 



L(q,Cl)6(u(u + Cl)) 



P 7 = 4ir 2 N(Q) 2 Tj2^2 



sgn(w + Q) 



(|w| + 1/t s )(ZV + |w| + |co + ft| + 2/t s )6(uj(uj + fi))) 



tfc(g,fi)L(g,n). 



(52) 



The assignment of terms to the polarization bubble diagram they would arise from if we had done the calculation by 
that method is unique, and can be summarised below: 

Pi : term proportional to Vc, 6(—lu(lj + Cl)) with no \u> + Cl\ + l/r s denominator. 

P2 : term proportional to Vc, 0(lj(uj + Cl)) with no \oj + Q| + l/r s denominator. 

P3 : term proportional to Vc, 9(—lo(lu + Cl)) with \u> + 0| + l/r s denominator. 

P 4 : term proportional to Vc, 0(lh(uj + Cl)) with no |u> + f2| + l/r s denominator. 

P5 : term proportional to L and o>(cj + Cl)). 

P 6 : term proportional to L and (9(o;(cj + CI)). 

P7 : term proportional to LVc- 
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We find that these reduce to the results of I when we set l/r s — > 0, providing a useful consistency check on the 
present calculation. 

To evaluate the correction to T c we split it into two parts: the Coulomb part consisting of those parts that contain 
a Coulomb propagator, (Pi — P4, P7), and consequently require special attention at q = 0, and the fluctuation part 
consisting of those terms that contain only a fluctuation propagator, (P5, Pg). Performing the w-sum first we get for 
the Coulomb part 



^(fe)- r -??{5P 

2Dq 2 p 2 + (Dq 2 ) 2 } 



Dq 2 



2ttT c n 2 - (Dq 



1 



1 



101 



2 2ttT c t s 



2TTT r 



|op 2 - (Dq 2 ) 2 } 2 



|0| (|o| + - 



|0| + Dq 2 



1 



101 



2 2irT c r s 2ttT c 



IMS 



1 



2ttT c t s 



4(Dq 2 ) 2 



101 



27rT r 



1 



2*kT c t & 



[Q2 _ (D g 2)2]2 



1 



2-kT c t r 



^g 2 + |0|) 



1 



1 



27rT n r« 



(53) 



Since the worst singularity possible in Vc(q, 0) at = goes as 1/q 2 , the overall g 2 factor multiplying Vc in the 
above expression means that this singularity is removed. It follows that the removal of the q — singularity in the 
Coulomb part is unaffected by the addition of magnetic impurities - this is because it is a general feature enforced 
by gauge invariance, as we will show in Appendix B. 

To calculate the Coulombic suppression term of Eqn. ([33]), we change variables to to = fl/2irT and y = Dq 2 /2irT, 
noting that V = J dy/(8n 2 DT), and 



N(Q)V c (q,n) = N(0) 

This leads to the result 
1 



2N(0)Dq 2 



Aire 2 Dq 2 



101 



Dq 2 + \Q\ m + y 



2Dq 2 



2?y 



(54) 



lo s( 



T, 



rO 



87T 2 iV(0)P> 



M 

E 



A/ 



(/(/ 



m,— y \Z 
2y(m 2 + y 2 ) 



2ay 



m(m 2 — y 2 ) 2 m(m + 2a)(m + y) 



ip [ — + at + m \ — i]j 



4y 



[^(\ + a + m) - (\ + a)\ 



(m~y)(m 2 ~y 2 ) ( I + a + 2±*) - (§ + a)] j ' 



(55) 



where a = 1/27tT c t s and the upper cutoff M — 1/2ttT c t. The leading order term is that which goes like l/y at large 
y, leading to logarithmic behavior. To isolate this, we add and subtract the term 



M 

E 



M 



dy 



2a 



' J (to + y) \ \ m m(m + 2a) 



1 



ip I — h a + to — t/j [ — ha 



1 



- , ' ( ~ + a 



M 

-E ln 

m— 1 

to give the result 
T c o 



M + m\ ( 2(m + a) 



m(m + 2a) 



+ a + to 



(56) 



,1/ 



Po > Z ^ 

M 

E 



m— 1 
M 



M 



2(m + a) 
m(m + 2a) 



■0 [ — + a + to 



_1 Jo 



m— 1 



4* 



+ a + m 



(to — y)(m 2 — y 2 ) 
1 



[ip(± + a + m) -ip(\+u)\ 



y-m (1 



(57) 
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where we have noted that l/8ir 2 N(0)D = Ru/Rq. We could now proceed to evaluate this expression, but before we 
do so, let us consider the domain of validity of the first-order perturbative result. 



V. BEYOND PERTURBATION THEORY 



Since we now have the full first-order perturbative correction to the transition temperature due to the effect of 
disorder on the Coulomb interaction, we could in principle plot curves of T c (R a , 1/t s ) and compare to experiment. 
However the curves of T c vs Ra for different values of l/r s would simply be exponential decays with different initial 
slopes. First-order perturbation theory is unable to treat the strong disorder region, and so cannot lead to the 
complete destruction of the superconductivity by non-magnetic disorder. If we are to consider the effects of arbitrary 
disorder strength, we must work beyond perturbation theory. In what follows we discuss two methods of doing this, 
and compare the results we obtain from them. 

The simplest way to proceed is to "self-consistently" solve the first-order perturbative expression of Eq. (|5^). 
This simply means that we replace T c q by T c on the right-hand side of Eq. (|5^), and solve the implicit equation we 
obtain for T r which has the form 



' [1 + 



1 



_^/(T c ,l/r s ). 



(58) 



Here f{T c , 1/t s ) is the complicated expression on the right-hand side of Eq. (57), whilst the first term is just the 
mean-field suppression of T c by the magnetic impurities. From our knowledge of the situation without magnetic 
impurities, we know that unphysical re-entrance problems may arise with the solution of this equation, and we should 
not take it too seriously in the region where superconductivity is strongly suppressed. 

The fact that we cannot trust the results obtained from this "self-consistent" theory leads us to ask the question 
of how to correctly go beyond first-order perturbation theory. The best approach is to deriveuthe effective field theory 
from which the perturbation series may be deduced - in this case a non-linear sigma modelE3 - and treat this using 
the renormalization group. This has been done by Finkel'steinEj for the system without magnetic impurities, but has 
the problem that it is_very difficult, and would become even more so if magnetic impurities were added. Recently 
Oreg and Finkel'steirO demonstrated that the same results could be obtained using a much simpler non-perturbative 
resummation technique, which we show diagrammatically in Fig. (7). The method uses a featureless Coulomb 
interaction of magnitude N(Q)Vc = 1/2, consistent with the cancellation of the l/q d ^ 1 divergence discussed earlier, 
and keeps only diagrams 3 and 4 of Fig. (6), since they give the greatest contribution. This leads to the equation for 
the pair scattering amplitude, T(Lu n ,uji), 



r(w n ,w/) = -|A| + tA(uj n ,uJi) - irT 



M 

E [ 

m=— (M+l) 



-|A| + tA(w n ,w m )]-. 



l/n 



-T(u! m ,U>l), 



(59) 



where u> n = 2-kT(ji + 1/2) is a Fermi Matsubara frequency, and the upper cut-off M — 1/2ttTt. The amplitude 
A(w„, lui) is given by 



A(uJ n ,LUl) = 



In 



hi 



1 



(\u) n \ + M)T 
1 



(\u n \ + H +2/r s )r 



LU n UJl < 
U) n UJl > 



(60) 



where the breaking of time-reversal invariance by the spin-flip scattering means that A has a different form depending 
upon the relative signs of its two Matsubara frequencies. If we treat the r(w„,w m ) as elements of a matrix T, the 
matrix equation for T can be solved to yield 



where 



f = Cj x i 2 (i - lAinr^-^HAii + tA), 



n = ^- 1 / 2 [i-|A|- 1 tA]^- 1 / 2 , 



(61) 



(62) 



tinm = (n + 1/2 + a)S nm , A nm = A(w n ,w m ), l nm = 1 and = 6 nm . The matrix T becomes singular when 
an eigenvalue of n equals 1/|A|, and this signals the onset of superconductivity. Note that the matrix II depends 
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on temperature both through the temperature dependence of its elements, and also through its rank 2M. To find 
T c , we start at the BCS value T c o, which corresponds to a value of M given by Mq — 1/2ttT c qt. We decrease the 
temperature T by increasing the upper cut-off M successively by one. For each value of M, we construct the matrix 
II, and diagonalise it. When its lowest eigenvalue equals 1/|A|, we have reached the transition temperature T c , which 
is given by T c /T c q — Mq/M. This method allows us to go to as low a temperature as we like, provided that we are 
prepared to diagonalize large enough matrices. 

We will now plot curves of T c vs Ra for fixed l/r s , and T c vs l/r s for fixed Ra, derived both from the self-consistent 
perturbation theory of Eqn. (|58|), and from the non-perturbative resummation approach of Eqn. (^9|). This is done in 
Fig. (8), and we see that the two approaches are in rough agreement. The resummation technique is seen to remove 
the re-entrance problem which occurs in the a = curve at large Ra, but surprisingly this re-entrance seems to be 
partially cured by the presence of magnetic impurities. 

The above curves are fine from the theorist's point of view, but experimentally what is measured is the suppression 
of T c by a certain fixed amount of magnetic impurities as the thickness of the superconductor is altered. The data is 
then presented in the form of the pair-breaking per magnetic impurity which can be written as 

a {Ra) = — , (63) 

which we can also generate from our theoretical expressions. The result will of course depend upon the magnitude of 
the value of a we choose: we would like to choose a as small as possible so that we are always in the linear regime 
of pair-breaking, but not too small so that the difference is very sensitive to the discrete sums used in the numerical 
calculation. A typical plot is shown in Fig. (9). If we ignore the numerical noise we see that a' is roughly constant and 
equal to its mean field value of 7r 2 /2. It only appears to increase as we approach the region where superconductivity 
is destroyed, and its total variation is only.-about 10% even if we include this region. This is in agreement with the 
experimental data of Chervenak and VallestJ. 



VI. DISCUSSION AND CONCLUSIONS 



The main conclusion of this paper is that the effects of localization and interaction do not lead to an appreciable 
change of pair-breaking rate per magnetic impurity in disordered superconducting films provided that we are not 
too close to the superconductor-insulator transition. The experimental data agrees with this theoretical prediction, 
and thus confirms the validity of the basic model of T c suppression in disordered superconductors which consists of 
the BCS interaction, Coulomb repulsion and static disorder. The fact that the theoretical prediction is obtained 
both from first-order perturbation theory, and from a non-perturbative resummation technique, gives us increased 
confidence in its validity. Our calculations demonstrate that the resummation technique is a very powerful tool for 
going beyond perturbation theory which can be adapted to a variety of situations. Moreover we find that the ad hoc 
"self-consistent" extension of first-order perturbation theory can give sensible results even at values of Ra near the 
superconductor-insulator transition, at least in the presence of pair-breaking. 

The effect of nonmagnetic disorder on pair-breaking in superconducting films has previously been considered 
by Devereaux and BelitzEll using a model in which strong coupling effects are considered. Good agreement with 
experiment is also obtained with this approach, although more fitting parameters are required in this model. We note 
that only a single fitting parameter - the initial slope of the T c (Ra) curve - is required in our approach. Unfortunately 
we see that the experimental data is unable to determine which, if any, of the two approaches is correct. In support 
of our approach we note that it is a "minimal model" in the sense that it contains the minimal physics to describe the 
system, and requires the input of a single fitting parameter. However, this is not to say that strong-coupling effects 
are not important in this system. 

Another important result which emerges from the approach based on the grand-canonical potential is that the 
1/q 2 singularity of the disorder screened Coulomb potential is always cancelled in first-order perturbation theory. 
This removes the possibility of changing some experimental parameter to obtain a strong suppression of T c from 
this singularity. We have shown that this cancellation is enforced by gauge invariance, and leads us to suspect that 
it occurs to all orders in perturbation theory. It is this cancellation which makes it legitimate to use a featureless 
interaction in the resummation technique. 

ACKNOWLEDGEMENTS 

We thank I. Aleiner, A.M. Finkel'stein and Y. Oreg for helpful discussions. R.A.S. acknowledges support from 
the Nuffield Foundation. V.A. is supported by the U.S. National Science Foundation under grant DMR-9805613. 



12 



APPENDIX A: CALCULATION OF POLARIZATION BUBBLES 



In this appendix we give a detailed derivation of the polarization bubbles, ILij, shown in Fig. (3). To evaluate 
these we must first calculate the impurity ladder, II, which is given by the geometric series 



where Tq is the impurity line 

r„= 1 



2ttN(0)t 



n = s + sr s + stqSTqS - 



[AiT 3 t7 ® T 3 f7 + A 2 (r cri (g) r ai + T a 2 ® t ct 2 + T 3 cr 3 ® r 3 cr 3 )] 



(Al) 



(A2) 



and 1/r = 1/to + l/r s is the total impurity scattering rate, Ai = t/tq, and A2 = t/3t s . S is the momentum sum of 
a direct product of Green functions 

S = ^2 G ( k > iuJ ) ® G ( k + ?) iuJ + 

(z- At 1( t 3 )(z' - A't iC t 3 ) 



7rW(0)r/ 



and / is the integral 



ee' 



6(6fe - q-v F ) 



(^-e 2 )[fe-q.v F )2-e' 2 ]' 



(A3) 



(A4) 



Since we do not need the complete matrix structure of II, but just its traces with two matrices from the set ti<t 3 , 
7203, 73CO1 we actually evaluate the impurity dressed vertices Tlj which have one matrix from the above set inserted 
between the two terms of the direct product in II. These satisfy the equation 



ILi = Si + STqIL 



(A5) 



Starting with IIa we see that 

5a = 27tAT(0)t/ 



T 3 (7oTlCT 3 T 3 £ro 



(z + ATl(T 3 )Tl(T 3 (z' + A TlCT 3 ) 



ee' 



-1 - 



= 2ttN(0)tI 

= 27riV(0)rJ(a + - /3 + tict 3 )ti(t 3 , 
where the a and (3 terms are coherence factors 

zz' ± AA' 



(z + Ana 3 )(z' + A na 3 ) 



ee' 



Tier 3 



a± = 1 



ee' 



a± = a± - 2 ; f3± 



z'A ± zA 



ee' 



By inspection we see that IIa must have the matrix form 

n A = 2ttN(0)tI[A + Bricr 3 ]riCT 3 , 



(A6) 



(A7) 



(A8) 



and we now substitute this into Eqn. (A5) to deduce the coefficients A and B. To derive the second term on the RHS 
of Eqn. (|a|) we see that 



ToIIa = —I\i(A - STicr 3 )ricr 3 + 3/A 2 (A + Sricr 3 )ricr 3 , 



(A9) 



and thus 
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5r n A = 2nN(0)rI L(A + Br^r^ + ^ ± + *™)™(* + 2™) 



££' 



-3A 2 (A - Br 1 (T3)ricr3 



(2 + Atio- 3 )\ 2 {A - Br 1 a z )T 1 az{z' + A na 3 ) 



= 2ttN{0)tI { (Ai - 3A 2 ) A 



(7+ An <T 3 ) (2' + A rio- 3 ) 



££' 



££' 



T1CT3 



-(Ai + 3A 2 )Sri(T3 



(z + Ana^iz' + A Tl a 3 ) 



££' 



no's 



27rjV(0)r/[(Ai - 3A 2 )A(-a + + p+na^r^ + (Ai + 3A 2 )£(a+ - /3 + Ti<7 3 )]. 



(A10) 



We can now equate the coefficients of 1 and T\a 3 on the LHS and RHS of Eqn. ( |A5| ) to obtain the linear equations 
for A and B, 



(All) 



This matrix equation can then be inverted by inverting the 2x2 matrix and using the identity a+a + = [3\ to obt am 



l + /(Ai 


- 3A 2 )a + 


J(A 1 + 3A 2 )/3+ 




A~ 




a+ 


-/(Ai - 


3A 2 )/3+ 


1 -7(Ai +3A 2 )a+_ 




B 




-0+ 



'A' 


1 


1 


B 


" ^7 







a+ 


1 


a+ 






~ ^7 





7(Ai-3A 2 )/3+ 1 + I(Ai - 3A 2 )S+ 

where the determinant D + can be written 

D+ = [l + /(Ax - 3A 2 )a+][l - Z(Ai + 3A 2 )a+] + I 2 (X 1 + 3A 2 )(Ai - 3A 2 )/?2 

+ AA /N 



= 1 + 7(Ai - 3A 2 )5+ - J(Ai + 3A 2 )a+ = 1 - 2/Aj + 6JA 2 



Similar results are obtained for II^ and n p , leading to the results 



££' 



(A12) 



(A13) 



n A - 2itN{0)t— (a+ - /5+T 1 a 3 )T 1( 73 
n = 2TrN{0)T^-(a- - 0-Tia 3 )T2a 3 
n p = 2ttN{0)t— (a_ - /3_t 1( 73)t3(7o, 



where 



D± = l- 2/Ai + 6/A 2 



z'± AA 



££' 



(A14) 



(A15) 



If we evaluate the integral in Eqn. (A4) we find that / is given by 

1 



2It = 



9 9 

q vf 



W + W 2(W + W 



where 



-7-2 



-i-'2 



W=\/lu 2 + A^ ; W=\Juj' 2 + A 
From the second part of Eqn. (|19|) we see that we can write 

W + W' = — - — + AU + AU', 

T T s 



(A16) 



(A17) 



(A18) 
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and substituting Eqns. (A18) and (A16) into Eqn. (A15) gives the result 
1 1 fuu' T 1' 



D+ = 1 



TO T S V UU' 



— - — + AU + AU' 

T T s 



Dq T : 



Dq 2 + AU + AU' 



, 1 / uu' =F 1 



T S V UU' 



T. 



(A19) 



We can finally obtain the non-zero polarization bubbles Ily by inserting the second matrix from the set T1CT3, T2&3, 
T3<7o into Tlj and taking the trace. This yields 



UJ 

u^{q : n) = ttN{o)tJ2 

n pp (q,n) = -irN(0)Tj2 

to 

U 4>p (q,n) = -7rN(0)TY^ 



UU' + uu'-l 



UU' 
UU' + uu' + 1 



1 



Dq 2 + AU + AU' 

T, 



UU' -uu' + 1 
UU' 



UU' 
UU' - uu' - 1 



Dq 2 + AU + AU' 



UU - uu' -1 
UU' 



UU' 



Dq 2 + AU + AU' 



UU' -uu' -1 



UU' 



N(0) 



UU' 



Dq 2 + AU + AU' 



UU' -uu' - 1 
UU' 



= -H^(q,Q). (A20) 



APPENDIX B: LOW-MOMENTUM SINGULARITIES IN DENSITY AND PHASE PROPAGATORS 



The identities 

- A" 1 + 11^(0, fi) = xU^ p (0, ft) , IT 0p (O, Q) - -xU pp (0, ft), (Bl) 
where x = f2/2A, play a central role in the present paper and in I. As we have seen, their consequence 

[-a- 1 + n 00 (o, o)]n w (o, q) + n 0p (o, nf = o (B2) 

leads to the potentials V^0, V^ p , and V pp having l/q^ 1 singularities at low momentum q for all temperatures < T < 
T c and all non-zero frequencies Q 7^ 0. The importance of these identities suggests that they embody an underlying 
invariance principle. In this appendix we show that they are Ward identities connected to charge conservation, which 
is very reasonable since the impossibility of instantaneously moving the conserved screening charge a finite distanea 
is at the root of these singularites, The ideas at work here go back to Nambu's 1960 paper and its elaborationsE3~E-3 
and they are only included here for completeness. It is unfortunate that the physical basis for these identities was left 
obscure in I. 

To avoid irrelevant notational complications, we shall work within the 2x2 Nambu space. The 4x4 space needed 
to deal with spin flip scattering does not affect the general argument, and we shall in any case explicitly verify the 
identities for this case later in this appendix. 

The 'proper polarization parts', II, are calculated within a mean field approximation, in which the interactions 
are replaced according to 

V -> Vmf = ATr^Vi*, (B3) 

which implies a choice of phase for the order parameter. It is known that the quasiparticles obtained in this ap- 
proximation do not conserve charge, because Vmf does not commute with the electron density. Since the only 
other non-commuting part of the Hamiltonian is the kinetic energy, the operator equation of motion for the density 
p = rj"*V 3 * is 

% + V • 3 = i[V MF ,p\ - 2ATWV 2 ir ; (B4) 
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where j = T? , [^V>I' — 'I'VvE'^] is the current density operator. Eq. (B4) leads to the identity 



S(xi - X2)5{h - t 2 )i[T 3 G(x2,t 2 ,x 3 ,t 3 )] lj - S(x 2 ~ x 3 )S(t 2 - t 3 )i[G(xi,t-L,x 2 ,t 2 )T 3 ] 

-r(*i(ii,ti)v 2 ■j(x 2) t 2 )<aUx 3 ,t 3 )) + T(* i (x 1 ,t 1 )p4x 2 ,t 2 )*Ux3,t 3 )), 



(B5) 



where T is the time ordering operator and we have defined p$ = Tr^T 2 ^ . [The first two terms on the right of Eq. 
(B5) come from the derivative of the time ordering operator.] Multiplying Eq. ( ]B5| ) by t 3 , taking the trace, and 
Fourier transforming in space and time leads in the limit of zero wave vector to the identity 



ttn pp (o,r!) = 2A n^(o,n). 

On the other hand, multiplying by t 2 and performing these same operations yields 



(B6) 



niWo.n) 



Tr [n G] + 2 A U 4 



= 2A[-A" 



if 



(B7) 



In th e se cond line above the self consistency equation within the mean field approximation has been used. To obtain 
Eq. (Bl) we must note that 11^ is antisymmetric in its indices because of time reversal invariance — under which p 
is symmetric and antisymmetric. 

Since the -impurity interaction commutes with the charge density these identities survive in any 'conserving' 
approximationcJ and, in particular, in our sum of all non-overlapping graphs. In the main body of this paper, the 
mean field approximation was used as a way station on the road to the single loop approximation of Section III. There 
the phase of ihe-prder parameter is not fixed as in Eq. ( B3 ) but determined self consistently, which restores charge 
conservation.c3n3 

Finally, we shall verify explicitly that the identities ( |Bl| ) are satisfied by our calculated expressions. We start 
with the equations for — A -1 + II^, Tl^p, and n pp , 

UU' + uu' + 1 



a- 1 +n^(o,n) = tJ2 
n 4ip (o,n) = TY / 



i 



UU' (U + U' -<i[ uu 'uw'~ 1 ]) u 
u' — u 

uu' (u + u' -c[ uu '-^'- i : \) 



n pp (o,ft) = --T^ 

7T * * 



UU' -uu' -1 



uu 1 (u + u> - C[ mi i%p L = 1 ]y 



(B8) 



where we have removed the common factor 7riV(0) from each II, and set A = 1 for algebraic convenience. These 
factors can, of course, be replaced when we have finished. 

We will first prove the relationship between — A -1 + IL^ and n^p, namely 



A- 1 + n^(o,n) = |$ 0p (o,o). 



To proceed note that we can write 

, 1, 



uu — — u 

2 L 



u' 2 - {u' - u) 2 ] = ^[U 2 + U' 2 - 2 - (u' - u)% 



(B9) 



(BIO) 



from which it follows that 



UU' +uu' + l = -[(U + U') 2 - (u' - u) 2 } 
UU' -uu' -1 = -i[(Z7' - U) 2 - (u' - u) 2 ]. 



In the last term on the RHS of Eq. (B8) for H^, we can use the transformation ui 
uu is invariant. This leads to«<-> — u and U «-> U' , so that 



(Bll) 

-uj', under which the sum over 
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2T y - = rV [- + — 1 = rV 



U + U' 
UU' ' 



(B12) 



We can now write, using Eqs. (B8), ( Bll ) and ( B12 ), 

'{U + U'f- {u' - u) 2 + VL{u' - u) U + U' 



2{-\~ l + n 00 ) - nn^ p = tJ2 



'UU'-uu'-l 



]) 



C/E/' 



[TCP (;/ + 17' - C L r^- 

(U + U'f - {v! - u)(u' - u- n) - (U + U') [U + U' - c 



UU'-uu'-l 



UU' 



uu' {u + u> -([ w-^'- 1 ]) 

From the definition of u and u' in Eqn. (|2C|) we obtain the identity 



(B13) 
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~ U 


» i 
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- + u- 


IF 



u u 

Tr^u 



u 



, 1 1 uu' uu' 
U + U 

u u> u w 



au + u r ) 



UU' -uu' -1 
UU' 



(B14) 



It follows that the numerator in Eqn. (B13) is zero, and hence we have proved the required result (B9). 
We next prove the relation between H pp and Ii^, p , namely 



n 0p (o,ft) = -^n pp (o,f>). 



We start by considering the sum 



lim T y 

un — >oc l- — J 



uj — — {uq+VL) 



a 



(B15) 



(B16) 



As |w| — > oo, u — > (\oj\ + C) s g n ( w ) an( l U /U — ► sgn(w). Since there are Q>/2ttT more negative terms than positive, the 
sum becomes 



lim T V 



u 

w 



V 2nT J V fi 



r 



i 

'2tt' 



(B17) 



cj — — (o;o + n) 

We have chosen the limits so that we are able to make the usual uj <-> — uj' transformation in this sum. It follows that 



1 = - 2T YJL = tY±-(?L-!)=tY 

m jjq O \TT> TT I Z^ 



1 /V 

n\u' u 



U + U' ( UU'-uu'-l 



n(u' - u) 



UU' 



(B18) 



where we first make the uj <-+ —uj' transformation, and then use (B14). We can then rewrite Eq. ( |BS| ) for H pp in the 
form 



k pp = tY 



UU' -uu' -I 
UU 1 



1 



U + U' 

n(u' - u )~ (U + U'-C [- V'uu" 1 ' 1 ] ) 



Ty, /UU'-uu'-l 



(u + u 1 ) [u + u' - c 



UU' 
UU'-uu'-l 



UU' 



(u-u')n 



(u'-u) (U + U'-C[ uu 'uu''^ ]) 



(B19) 



From the identity (B14) we see that the numerator of (B19) can be rewritten as 
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(U + U') [u + u'-c 



UU' -uu'-l 



uw 



J - (V - u)fi 



= (c/ + c/') 2 - (V - «)(«' - u - fi) - K - «)fi 

= (C/ + C/ , ) 2 "(«'"«) 2 , 



and inserting the identity (|B14| ) into the first factor in ( B19 ), we get 



n 



pp 



2VL 



E 



(u> - u)UW (U + U'-C [ ""'uw'" 1 ] ) 



Multiplying out the numerator yields 



from which it follows that 

n„(o,n) 



[{U 1 - Uf - (u' - u) 2 ][(U + U'f - [v! - uf] 
(u' 2 - u 2 ) 2 - 2(u' - uf(U 2 + U' 2 ) + («' - uf 
{v! - uf[(u + u'f - 2u 2 - 2u' 2 - 4 + («' - uf] 
-i(u'-uf, 



^ --in (no) 

' n ^ uu' (u + u' - t fJV'-^'- 1 }) n M ' J ' 



(B20) 



(B21) 



(B22) 



(B23) 



completing our proof of the result (B15) 



APPENDIX C: EVALUATING DERIVATIVES OF Uu WITH RESPECT TO A 



In this appendix we evaluate the derivatives of the polarization bubbles Ily with respect to the order parameter 
A so that we may evaluate th e firs t order correction to the order parameter self-consistency equation. The formulas 
for the Hij are given in Eqn. ( A20 ), and we see that the derivative can operate either on the coherence factor or the 
denominator present in these expressions. 

The difficulty in evaluating these derivatives arises because u(uj) satisfies the transcendental equation 



A 



1 



1 



At, (u 2 + 1)V2 



(CI) 



from which it follows that 



A 2 " 



du 
OA 



1 - 



Ar s (u 2 + l) 3 / 2 



A 2 t s (t 



1)1/2 ' 



(C2) 



and thus 



du 
OA 



u 
A 



Ar s (u 2 + l) 3 / 2 



(C3) 



with a similar result for v! . 

We first consider the effect of d/dA on the coherence factors present in the Hy. We see that it suffices to evaluate 
the derivatives 



8 
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u — u 
UU' 



(C4) 



The first term in Eqn. (C4) gives 
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u 3 



u'). 



(C5) 
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Since this expression will occur inside a sum over u, and will multiply an expression that is invariant under the 
transformation lu <-» —u>', we see that the two terms in Eqn. (C5) will give equal results. Thus 



d_ 
dA 



UU' 



2 uu' 
AU 3 U' 



u 3 



The second term in Eqn. ( |C4|) gives 

1 



_d_ 

dA 



UU' 



1 u du . ,. 2 u 2 

-wTp^ + {u ^ u) = auW 



whilst the third term in Eqn. (C4) gives 



_d_ 
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UU' 



u' u du 1 1 du , _ 2 u(uu' + 1) 

TriPdA ~ TriPdA + u ' = a u 3 U' 



! -1 



1 - 



u 3 



(C6) 



(C7) 



(C8) 



The effect of d/dA on the coherence factors can then be summarised in the form 



d_ 
dA 
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dA 



uu' ± ll 2 u(u' =F u) 

UU' \ = _ A U 3 U' 
u' — u 2 u(uu' + 1) 

UU' = A U 3 U' 



U 3 



U 3 



Next we must consider the effect of d/dA on the two denominators 



D± = Dq 2 + AU + AU' 



uu' =F 1 
UU' 



(C9) 



(CIO) 



The last term on the RHS is a coherence factor, and its derivative can be read off from Eqn. ( p9|) above. The only 
terms then to consider are AU and AU' , which, of course, will give identical results after summation over u>. We see 
that 



d , . TTX TT . u du u 2 
M (m = U + A^- = U-- 3 
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From this we obtain the final result 
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U U 2 
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(Cll) 



(C12) 



Having now evaluated the action of d/dA on all the components of the polarization bubbles, Ely-, we can now 
write down the results for the dUij/dA, 
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(C13) 
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FIGURES 




X 



(a) (b) 

FIG. 1. The electron Green function in a superconductor, (a) Self-energy for a clean superconductor. The wiggly line is the 
BCS interaction, (b) Extra self-energy diagram needed for dirty superconductor. 
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FIG. 2. Definition of the screened potential Vij in terms of the polarization bubble 11^ and bare potential V°j. 
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(a) 




FIG. 3. Definition of the polarization bubbles IT^. (a) The geometric series for the ladder II. (b) the geometric series for the 
vertex function Hj which is obtained from IT by taking the trace at one end with a Pauli matrix, (c) The polarization bubble 
Hij is obtained from the vertex operator IT-, by taking the trace with a Pauli matrix at the open end. 
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FIG. 4. The first-order correction to the grand canonical potential. This has the form of a "string of bubbles", where the 
wiggly lines can be either the bare Coulomb or BCS interaction, and the bubbles are any of the non-zero polarization bubbles. 
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(a) (b) (c) 

FIG. 5. The equivalent Eliashb erg- like self-energy diagrams for the correction to the order parameter A. (a) The two-ladder 
diagrams are obtained by differentiating the diffusion propagator term in the Tlij with respect to A. (b) The one-ladder 
diagrams are obtained by differentiating the coherence factor term in IT-, with respect to A. (c) No three-ladder terms are 
obtained by differentiating fii(A) with respect to A, and direct calculation of these diagrams shows that they equal zero. 
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(a) 




(la) (lb) (lc) (2) 




(4b) (5a) (5b) (5c) 




(6) (7) 

(c) 

FIG. 6. The first-order correction to the pair propagator. The value of temperature at which this first diverges is the 
transition temperature T c . (a) Definition of pair propagator in terms of BCS interaction and pair polarization bubble, (b) 
Zeroth-order (mean field) pair polarization bubble, (c) The 7 diagrams which contribute to the first-order correction to the 
pair polarization bubble. The wiggly line is the screened Coulomb interaction, whilst the spring-like line is the pair propagator, 
as defined in (a). 
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FIG. 7. Diagrammatic equation for the scattering amplitude matrix T(uj„,lji). Block 7 is the BCS interaction. Block tA is 
the correction to the effective interaction caused by the interplay of Coulomb interaction and disorder. 
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Resistance R/£l aJa 

FIG. 8. Plots of transition temperature as a function of resistance per square and spin-flip scattering rate. The plot on the 
left shows T c as a function of Bu for values of (top to bottom) a = 1/2tvT co t s equal to 0, 0.2, 0.4, 0.6 and 0.8 times the critical 
value ao. We see that the a — curve has a re-entrance problem, but that the situation improves for finite a. The circles are 
the results from the non-perturbative resummation technique. We see that they are roughly in agreement with the perturbation 
theory. The plot on the right is of T c as a function of spin-flip scattering measured in the dimensionless form a/ao for values 
of Rn equal to (top to bottom) 0f2, 500J1, 1000fi, 1500O and 2000fX The circles are the non-perturbative resummation results, 
and are again in good agreement with perturbation theory except for the 2000fi curve. This might be expected since this curve 
is very close to the superconductor-insulator transition. 
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FIG. 9. Plot of pair-breaking rate per impurity versus resistance per square of film. We see that this is roughly constant, 
only increasing very near to the superconductor-insulator transition, with a variation of only 10% over the whole range. The 
curve on the left is from perturbation theory; the curve on the right from the non-perturbative resummation. 



29 



